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We show that the magnetic moment of a composite multiferroic torsional oscillator can be switched 
by the electric held. The 180° switching arises from the spin-rotation coupling and is not prohibited 
by the different symmetry of the magnetic moment and the electric held as in the case of a stationary 
magnet. Analytical equations describing the system have been derived and investigated numerically. 
Phase diagrams showing the range of parameters required for the switching have been obtained. 
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I. INTRODUCTION 

There has been enormous interest in investigating the 
switching of the magnetization of a nanomagnet by the 
spin-polarized electric current or by the electric held as 
possible routes to increasing computer speed. Spin trans¬ 
fer torque random access memory (STT-RAM) devices 
have been successfully engineered as cache memory and 
have entered the market. Large effort has been developed 
to achieve fast switching of the magnetic moment by the 
electric held in magnets interfaced with piezoelectric ma¬ 
terials and in multiferroic systems m ■ This latter task 
is difficult to achieve due the linear relation, M ocE, be¬ 
tween the magnetic moment M and the electric held E is 
prohibited by symmetry. It prompted researchers to look 
for, e.g., 90-degree rotation of the magnetization through 
the effect of the electric held on the magnetic anisotropy, 
or consider combined effects of the electric held and the 
exchange bias in layered systems. 

In this paper we propose a nanoelectromechanical sys¬ 
tem (NEMS) for switching of the magnetic moment by 
the electric held that is based upon a different principle. 
It follows from the observation that symmetry permits 
linear relation between M and the angular velocity 
of the rotation. The latter can be generated by a pulse 
of the electric held applied to a multiferroic particle at¬ 
tached to a torsional cantilever, or by the electric held 
pulse applied to a ferroelectric cantilever containing a 
magnetic particle. According to the Larmor theorem, in 
the coordinate frame of a nanomagnet, torsional vibra¬ 
tions are equivalent to the ac magnetic held 00. The 
latter is known to be capable of switching the magnetic 
moment, see, e.g., Ref. 0 and references therein. The 
advantage of mechanical vibrations over the ac magnetic 
held is that they can be localized at the nanoscale. This 
permits miniaturization down to molecular dimensions 
as has been demonstrated by recent measurements of a 
single molecular spin in a NEMS obtained by grafting of 
a single-molecule magnet on a carbon nanotube HOI HU- 

Coupling of mechanical vibrations to the magnetiza¬ 
tion of a thin him deposited on a small cantilever has 
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FIG. 1: Schematic representation of the geometry of the 
model studied in the paper. A particle with a magnetic mo¬ 
ment M and electric polarization P is a part of a torsional 
oscillator rotating about the X axis. The system is placed in¬ 
side a pair of capacitors. The reversal of the magnetic moment 
is generated by pulse voltages, V x and V y , on the capacitors. 


been studied before [T2] . In a multidomain magnetic him 
it is related to the motion of domain walls p~3l4T5j . In a 
single-domain system the magnetomechanical magnetiza¬ 
tion reversal was shown to have strong sensitivity to the 
parameters m ■ The concept of spin-rotation coupling 
sheds new light on the origin of the observed stochastic- 
ity. Circularly polarized ac held switches the magnetic 
moment in a predictable manner via the consecutive ab¬ 
sorption of photons of one helicity (i.e., one spin pro¬ 
jection) while a linearly polarized held consists of pho¬ 
tons of different helicity that are absorbed intermittently, 
raising or decreasing the corresponding projection of the 
magnetic moment. Although torsional vibrations cannot 
be reduced to either linear or circular polarized ac held, 
understanding of how the ac held reverses the magnetic 
moment [9: allows one to design the optimal geometry in 
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which rotations induced by the electric field can do the 
job. 

The equivalence of mechanical rotations to the ac mag¬ 
netic field in the rotating frame makes the problem of 
the electromechanical magnetization reversal uniquely 
defined and free of any unknown coupling constants. In¬ 
teraction of the magnetic moment with torsional oscil¬ 
lations is completely determined by the Larmor theorem 
and by the dynamics of the total angular momentum m 
spin plus mechanical, in the presence of the mechanical 
torque generated by the electric field. The paper is orga¬ 
nized as follows. The system under consideration and the 
model are described in Section |TI| where dynamical equa¬ 
tions for the magnetic moment and mechanical rotations 
caused by the time-dependent electric field are derived. 
Numerical solution of the equations and switching phase 


diagrams are presented in Section |III| Estimates for real 
systems and possible applications of our results are dis¬ 
cussed in Section El 


II. THE MODEL AND DYNAMICS OF THE 
SYSTEM 


No interaction of M with P in addition to the magnetic 
anisotropy is assumed. The effective coupling between 
M and E is provided by the rotation of the magnetic 
anisotropy axis that follows the direction of P. The eas¬ 
iest way to take this effect into account is by noticing 
that according to the Larmor theorem the rotation with 
the angular velocity ft = (dp/dt)x generates the effective 
magnetic field 


TT _ 1 d P~ 

7 at 


(3) 


in the coordinate frame of the rotating particle. In what 
follows we shall assume that no external magnetic field 
is applied. Correspondingly, H# is the only field felt by 
the particle in its own coordinate frame. 

While the easy magnetization axis in our model follows 
the direction of the electric polarization P, the latter does 
not necessarily follow the direction of the electric field E 
because of the finite moment of inertia of the rotator. We 
shall now proceed to the derivation of coupled equations 
for the magnetic moment and mechanical rotation caused 
by the application of the time-dependent electric field. 


A. Model 

We consider the system schematically represented in 
Fig. [TJ Magnetic energy of the uniaxial ferromagnetic 
particle in its own coordinate frame is given by 

E m = -KVm 2 z - Vm ■ H, ( 1 ) 

where K is the magnetic anisotropy constant, V is the 
volume of the particle, H is the magnetic field in the 
frame of the particle, and m = M/E is its constant- 
length magnetization. The body-frame Z-axis in this 
formula coincides with the easy magnetization direction. 
The particle is allowed to rotate mechanically around the 
X-axis, with p being the angle that the easy axis makes 
in the YZ-plane with the stationary Z direction in the 
laboratory frame. 

We shall assume that, in addition to the magnetic 
moment, the system also possesses the electric polariza¬ 
tion characterized by the vector P directed along the 
anisotropy axis of the particle. This can be achieved by 
either using a composite multiferroic material or by in¬ 
terfacing the magnet with the ferroelectric. Interaction 
of the electric polarization with the crystal lattice, being 
of Coulomb origin, is strong, providing that the orien¬ 
tation of the particle in space follows the orientation of 
P. On the contrary, interaction of M with the mechani¬ 
cal rotation of the particle is of the relativistic spin-orbit 
origin. Consequently, the magnetization vector does not 
necessarily follow the rotations of the anisotropy axis if 
the latter are sufficiently fast. 

Mechanical rotation of the particle in our model is 
achieved due to the interaction of P with the electric 
field E, 

( 2 ) 


B. Dynamical equation for the magnetic moment 
in the rotating frame 


The equation of motion for s = m/|m| in the coordi¬ 
nate frame rigidly coupled to the particle is m 

ds 

— = 7 [s x H eff ] - cry [s x [s x H eff ]], (4) 

where 7 is the gyromagnetic ratio, a is the dimensionless 
damping coefficient, and H e ff is the effective magnetic 
field given by 

HeS = ~V^ = HA + HR ’ (5) 


with H a = 2Km z z being the effective field due to the 
magnetic anisotropy. 

Introducing the frequency of the ferromagnetic res¬ 
onance, cjfmr = 2 yifm, and the dimensionless time 
t' = cjfmrC Eq. Q can be written as 


ds 

dt' 


s x 



dpjt^ A 

dt' ) 


a 


s x 


s X \s z z + 


dp(t') 

dt' 


x 


(6) 


At small damping, a 1, the last term in Eq. © can 
be replaced with — a[s x ds/dt']. 

Writing s in spherical coordinates in the body frame 
of the rotator, s = (sinO cos </>, sin 0 sin 0 , cos 0 ), one ob¬ 
tains the following equations for the angles describing the 
direction of the magnetization with respect to the fixed 
axes inside the rotating particle 


dO 

dt' 


dp . . . . dp 

— sin <fi — a sin 0 cos 0 + a —— cos 0 cos <fi. 
dt' dt' 


( 7 ) 


E P = -P E. 
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d(\> 

dt' 


„ dp 

cos 6 + — cos 0 cot 6 — a 
dt' 


dp sin 0 
dt' sin 0 ’ 


(8) 


Before studying rotations generated by the electric 
field, one can ask what p(t) would provide the optimum 
magnetization reversal. In the case of the ac field such 
optimum reversal occured in one plane, 0 = const 
Setting 0 = (po in Eqs. 0 and 0 one obtains at a = 0 

dp 1 dO sin# 
dt' sin po dt' cos po ’ 

Integration gives 


Here uo r is the resonance frequency of the torsional os¬ 
cillator. Note the proportionality of the elastic energy 
to p 2 . In nanowires it is justified even for large rotation 
angles m- The torque is given by 


t x 


dU 

dp 


— I r id 2 p + 


rf(P • E) 
dp 


(15) 


where P = P(0, — sin p, cos p). 

For, e.g., the electric field directed along the F-axis, 
E = (0, E, 0), one obtains from Eq. (12) the following 
dimensionless torque equation 


tan ^ = e*’ tan Vo ( 10 ) 

which corresponds to the switching of the magnetization 
from 0 = 0 at t' =» —oc to 0 = tt at t' = Too. According 
to Eq. 0 the function p(t') that generates such behavior 
is given by 


2 arctan[exp(t'tan (po)] 

^ sin <p 0 ' 

It describes the rotation that exponentially turns on at 
p = 0 and dies out at p = 7r/sin<po, thus, effectively a 
characteristic time interval At 1/(wfm,r tan</? 0 )- This 
kind of switching requires the rotation of the magnet by 
a significant angle. Curiously, if po = 30°, the 360° rota¬ 
tion does the job. It provides a good hint into the speed 
and amplitude of the rotations that should be imple¬ 
mented with the help of the electric field pulse to achieve 
magnetization reversal. 


C. Dynamical equation for the torsional oscillator 


Modeling of the effect of the field requires knowledge 
of the equation of motion of the torsional oscillator that 
must be solved together with the dynamical equations 
for the magnetization, Eqs. ([T]) and (| 8 |. For rotations 
about the X-axis it follows from the equation 


dJx 

~dt = Tx ’ 


( 12 ) 


where J is the total angular momentum and r is the 
mechanical torque. The total angular momentum con¬ 
sists of the spin angular momentum, hS = M /7 and the 
mechanical angular momentum, I r dp/dt , 

j x = ns x + i r ( 13 ) 

where I r is the moment of inertia of the rotator. 

When the rotation is generated by the electric field E 
acting on the electric polarization P, the potential energy 
of the rotator consists of the elastic part and the part 
generated by the interaction with the field, 


U = Ee + Ep 


1 

2 


i t < 4 p 2 - p • e. 


(14) 


[p+p{t’) COS P \ = 




where 


= 


UJ r 


UFMR 


II = 


Ir^FMR 

ns ’ 


pit') = 


PE(t') 

irW2 ' 


(16) 


(17) 


This case corresponds to the utilization of only the F- 
capacitor shown in Fig. [1] Notice that we have intro¬ 


duced the mechanical damping term in Eq. (16), propor¬ 
tional to the damping parameter 77 . 

For torsional oscillations about a fixed axis, induced by 
a pulse of the linearly polarized electric field, as in Eq. 
(16), the problem resembles the problem with a linearly 
polarized ac magnetic field for which the chaotic behav¬ 
ior of s should be expected. It makes sense, therefore, to 
consider a two component electric field, E = (0,£0,E0), 
generated by two pairs of capacitors shown in Fig. [l] 
One has to find the optimal E (t) that reverses the mag¬ 
netization with the initial conditions E = 0 , # = p = 0 , 
at t = 0. The most promising case turns out to be the 
case of a rotating field: E = E(t)( 0, — sin oo't', cos u't') of 
varying amplitude E(t). Then P E = PE(t ) cos(p-uj't') 
and the torque on the magnetic particle is 


-I r uJrP - 


PE(t) sin(p - 


ou't') 


(18) 


Instead of Eq. (16) one obtains 


dt^ + [P + P(0 sin (# “ dt 1 1 

= -7T5? (sin,W) 


(19) 


The right-hand-side of this equation represents the Ein¬ 
stein - de Haas effect, that is, the back effect of the change 
in the magnetization on the mechanical motion of the os¬ 
cillator. It is small when parameter representing the 
ratio of the mechanical angular momentum and the spin 
angular momentum, is large. Notice that a stationary 
electric field directed along the stationary Z-axis in the 
laboratory frame corresponds to u/ = 0. In this case 
p = 0 , s = const is the solution of the equations in ac¬ 
cordance with the expectation that such a field does not 
generate mechanical rotation. 
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FIG. 2: Switching dynamics of the magnetization in the co¬ 
ordinate frame of the magnet plotted together with its rota¬ 
tional dynamics for two sets of parameters. Blue line shows 
the projection of the normalized magnetic moment on the 
magnetic anisotropy axis. Yellow line shows spatial orienta¬ 
tion of the torsional oscillator given by p/(2tt). 


FIG. 3: a) Dynamics of the magnetization and rotational dy¬ 
namics of the oscillator showing no switching at long times for 
the chosen set of parameters, b) Switching of the magnetic 
moment accompanied by the permanent twist of the oscilla¬ 
tor (p ^ 0) at the end of the field cycle for another set of 
parameters with po > p c - 


III. NUMERICAL RESULTS 


In this Section we present results of the numerical 
study of the magnetization reversal governed by equa¬ 
tions ([7]), ([ 8 ]), and ( p~9] ) . In the numerical work we choose 
a = 0.01 for the Landau-Lifshits-Gilbert damping pa¬ 
rameter, which is in the right ballpark for single-domain 
magnetic particles [ 20 ]. As long as a is small, our results 
are not very sensitive to its exact value. They exhibit 
much greater sensitivity to the damping of the oscilla¬ 
tions of the torsional resonator p. Strong mechanical 
damping provides deterministic ranges of parameters for 
which the switching occurs, while weak damping gives 
more stochastic dependence on parameters. For that rea¬ 
son we present numerical results for two cases: Weak 
mechanical damping, 77 = 0 . 01 , and strong mechanical 
damping, 77 = 0.5. 

We apply a pulse of the rotating electric field, E = 
Eq ( 0 , — sinc^'t', coscc/t'), of constant amplitude Eq to the 
resonator shown in Fig. [lj The field is initially aligned 
with the electric polarization and then is rotated by 360°, 
which corresponds to a non-zero value of uj f within the 


time interval satisfying, 0 < uj't' < 2i r and zero uj' out¬ 
side that time interval. This generates mechanical oscil¬ 
lations of the resonator and the ac magnetic field in the 
coordinate frame of the magnet. From practical point of 
view one only needs to switch the magnetization, preserv¬ 
ing the orientation of the magnetization carrier in space. 
Static solutions of Eq. (19) satisfy 


p + po sin p = 0 


( 20 ) 


The return of the oscillator to its initial, p = 0, orienta¬ 
tion in space at the end of the field rotation is guaranteed 
if p = 0 is the only solution of Eq. (20). It is easy to see 
that this occurs when the amplitude of the electric field 
satisfies Po < p c = 4.60334. At high fields the oscillator 
stays twisted (that is, with p ^ 0 ) at the end of the field 
rotation if the field is not switched off or reduced below 
Po = Pc- 


Typical dynamics of the magnetization and the rota¬ 
tional dynamics of the oscillator are illustrated by Fig. 
[2] and Fig. [3] The dynamics is sensitive to the ratio, 
uj f rl of the resonant frequency of the oscillator and the 
FMR frequency, and to the ratio, u/, of the frequency 
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FIG. 4: Switching phase diagram in the phase space (oj',oo' r ). FIG. 5: Switching phase diagram in the phase space (po,uj) 
Areas of magnetization reversal are shown in blue (dark) representing the effect of the rotating electric field. Areas of 

color. magnetization reversal are shown in blue (dark) color. 


of the field rotation and the FMR frequency. It is also 
sensitive to the electromechanical factor po that is pro¬ 
portional to the amplitude of the electric field Eq, and to 
the mechanical damping of the oscillator 77 . Switching vs 
non-switching is less sensetive to the parameter a when it 
is small, and has little sensitivity to the normalized mo¬ 
ment of inertia of the oscillator I' r as long as it is large, 
meaning that under these conditions the damping term 
in Eq. 0 and the right-hand-side of Eq. ( |l9| ) contribute 
weakly to the dynamics of the system. In the numeri¬ 
cal work we choose I' r = 100. Magnetization reversal for 
strong and weak damping of the oscillator is illustrated 


in Fig. [2] For a different set of parameters chosen in Fig. 
it no reversal takes place at the end of the field-rotation 
cycle. When po > p c the oscillator, in accordance with 
expectation, ends up in a twisted state with a non-zero 
asymptotic value of p, see Fig. Ik 

f^. a shows the switching phase diagram in the 
(a/,a;') plane. Weak mechanical damping, 77 = 0.01, re¬ 
sults in a stochastic dependence on parameters, Fig. Ik 
which is a consequence of the high non-linearity of the 
dynamical equations. Strong damping, 77 = 0.5, gives 
more deterministic behavior characterized by a compact 
area of the magnetization reversal, Fig. it- Both types 
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of behavior exhibit thresholds on uo f and uj' r required for 
the reversal. 

Switching phase diagram in terms of the parameters 
representing the rotating electric field, (po , cU), is shown 
in Fig. [5] Here again the weak mechanical damping of 
the oscillator results in a stochastic dependence of the 
switching on the parameters, which is illustrated by Fig. 
UK while strong damping gives more deterministic behav¬ 
ior, Fig. If- It shows islands in the (po,cj') plane that 
provide the required strength of the electric field and the 
speed of the field rotation. 

IV. DISCUSSION 

We have shown that a magnetic moment of a single¬ 
domain magnet that is allowed to perform torsional os¬ 
cillations with respect to a fixed axis can be switched by 
the mechanical oscillations alone. This effect arises from 
the fact that in the reference frame of the magnet the 
rotation is equivalent to the magnetic field, thus mak¬ 
ing torsional oscillations equivalent to the ac magnetic 
field. It is well known [9] that the ac magnetic field, 
even of small amplitude, is capable of reversing the mag¬ 
netic moment. The physics of the effect is related to the 
consecuitive absorption of photons that, having a non¬ 
zero angular momentum, change the projection of the 
spin angular momentum of the magnet. In that sense 
the mechanical rotations are equivalent to photons, with 
the spin angular momentum of the magnet absorbing the 
mechanical angular momentum associated with torsional 
oscillations. 

In this paper we propose to ignite the required me¬ 
chanical oscillations by the time dependent electric field 
(voltage) applied to a multiferroic system or to a ferro¬ 
electric firmly coupled with a magnet. According to Fig. 
[4j the switching of the magnetic moment occurs when the 
frequency characterizing the rate of change of the field, 
and the resonant frequency of the oscillator, are not very 
small compared to the frequency of the ferromagnetic res¬ 
onance due to magnetic anisotropy (no external magnetic 


field is required). The latter is typically in the ballpark 
of a few GHz. Nano mechanical oscillators of such fre¬ 
quencies are common nowadays. (Oscillation frequencies, 
/ = uj/{2i r), of hundreds of GHz have been reported in 
NEMS based upon carbon nanotubes [21].) The effective 
magnetic field due to torsional oscillations of amplitude 
Po at a frequency uj = 2irf is of order h rsj upo/l- For / 
of a few GHz and po ~ 0.1 it is in the balpark of lOmT, 
which is more than sufficient for magnetization reversal 
by the ac magnetic field [9], thus allowing for even weaker 
torsional oscillations than po ~ 0.1 used in the estimate. 

The effect of the electric field on the oscillator is de¬ 
scribed by the electromechanical ratio po = PEo/(/ r co^). 
Fig- [5] shows that the field pulse results in the magne¬ 
tization reversal when p 0 is not too small compared to 
one. Estimating the electric polarization P as PoV and 
the moment of inertia of the oscillator as I r ~ mr 1 2 3 4 * , with 
Po, V, m and r being the polarization density, the vol¬ 
ume, the mass, and the size of the oscillator, respectively, 
one obtains that the voltage, Ve ~ Eor required to op¬ 
erate the device must be of order Ve ~ moj^/Po- For 
m ~ 10 ~ 20 kg, uj ~ 10 10 s -1 , and Po ~ 0.1 C/m 2 , this 
gives voltages in the ballpark of a few volt. Making such 
a device, while challenging, must be a doable task in the 
light of the recent progress in manufacturing NEMS. No¬ 
tice that high quality factor of the nanooscillator that is 
required for their applications as sensors is not needed for 
the purpose of magnetization switching. In fact, as our 
study shows, strong mechanical damping of the oscillator 
produces more reliable results. 
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